Disorder-assisted melting and the glass transition in amorphous solids 
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The mechanical response of sohds depends on temperature because the way atoms and molecules 
respond collectively to deformation is affected at various levels by thermal motion. This is a funda- 
mental problem of solid state science and plays a crucial role in metallurgy, aerospace engineering, 
energy. In glasses the vanishing of rigidity upon increasing temperature is the reverse process of the 
glass transition. It remains poorly understood due to the disorder leading to nontrivial (nonafiine) 
components in the atomic displacements. Our theory explains the basic mechanism of the melting 
transition of amorphous (disordered) solids in terms of the lattice energy lost to this nonaffine mo- 
tion, compared to which thermal vibrations turn out to play only a negligible role. It predicts the 
square-root vanishing of the shear modulus G ^ \/Tc — T at criticality observed in the most recent 
numerical simulation study. The theory is also in good agreement with classic data on melting of 
amorphous polymers (for which no alternative theory can be found in the literature) and offers new 
opportunities in materials science. 
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The phenomenon of the transition of a supercooled liq- 
uid into an amorphous solid has been studied extensively 
and many theories have been proposed in the past, start- 
ing with the Gibbs-DiMarzio theory fV . All these theo- 
ries focus on the fluid to solid aspect, coming to this glass 
transition from the liquid side (supercooling). However, 
it is only one facet of the problem. For the reverse pro- 
cess, i.e. the melting of the amorphous solid into a liquid, 
no established theories are available. 

The problem of describing the melting transition into a 
fluid state [2]-[9] is complicated in amorphous solids by the 
difficulties inherent in describing the elasticity down to 
the atomistic level (where the thermal fluctuations take 
place). It is well known that the standard (Born-Huang) 
lattice-dynamic theory of elastic constants, and also its 
later developments [TO^, breaks down on the microscopic 
scale. The reason is that its basic assumption, that the 
macroscopic deformation is affine and thus can be down- 
scaled to the atomistic level, does not hold [11]. Atomic 
displacements in amorphous solids are in fact strongly 
nonaffine pTHMi , a phenomenon illustrated in Fig. [l] 
Nonafiinity is caused by the lattice disorder: the forces 
transmitted to every atom by its bonded neighbors upon 
deformation do not balance, and the resulting non-zero 
force can only be equilibrated by an additional nonaffine 
displacement, which adds to the affine motion dictated 
by the macroscopic strain. 

Recently, it has been shown [TF that nonaffinity could 
play a role in the melting of model amorphous solids, al- 
though the basic interplay between nonaffinity, thermal 
expansion, and thermal vibrations remains unclear. With 
a number of other models of the glass transition, such as 
mode-coupling theories, there is also an issue when they 
rely on liquid-state theory for relating the stress tensor to 
local fluctuations in the solid state. These theories pre- 
dict that the shear modulus G of athermal hard-sphere 
colloids remains finite at the glass transition and 



that it jumps discontinuously to zero upon decreasing the 
packing fraction cj). This scenario does not agree, how- 
ever, with thermal systems, where both simulations [17] 
and experiments [T8^ show that the vanishing of G is 
continuous with T. The glass transition temperature Tc 
could therefore be estimated for thermal systems using 
the Born melting criterion [2], i.e. setting G(T) = [19j 
and solving for Tc, if a theory for the low- frequency G in 
the solid state is available. 

Here we discover a much more basic effect that non- 
affinity of deformation causes in the melting, where 
marginal stability controls the response through the 
bond-connectivity, in turn affected by both T (via ther- 
mal expansion) and disorder. The resulting framework 
allows us to link lattice disorder, connectivity and ther- 
mal vibrations into the first description of the classical ex- 
periment [20 on T-dependent low-frequency shear mod- 
ulus of amorphous polymers. 

Let us start from the basic assumptions of Born-Huang 
lattice dynamics [2 . The free energy density of affine de- 
formation is given by the following harmonic lattice sum: 
= w^ij (^^^/^^ij)i?o which runs over all 

bonded atom pairs ij. Here U{r^j) is the pair-interaction 
potential and the vector ufj = rfj — R^j = 77 • R^j de- 
notes the affine displacement, with 77 the macroscopic 

strain tensor and R,, the bond vector in the undeformed 

— ''J 

frame. Rq is the equilibrium lattice constant in the unde- 
formed frame, at which one evaluates the lattice spring 
constant k = (d'^U/dr^j)^ . Without loss of general- 
ity, we focus on shear strain r]xy = 7. The lattice sum 
can be evaluated upon introducing the average number of 
bonds per atom and in the afiine approximation [2T] : 
Fa = ^^{tv/Ro)(l)zj'^ . Here (j) = vN/V is the packing 
fraction occupied by the atoms or building blocks of the 
solid. 

In general, the contribution to F due 
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FIG. 1: a, If the deformation was affine, the particles labeled 
with darker color would still lie on the dashed lines also in 
the deformed (right) frame. Because of nonaffinity, they do 
not. b, In an ordered crystal lattice (left) the forces (arrows) 
transmitted to any particle i by its nearest neighbors upon 
deformation cancel each other and the net force acting on i 
is zero by symmetry. In a disordered lattice (right) the forces 
transmitted upon deformation do not balance, and a net force 
acts on i which has to be balanced by a nonaffine motion, to 
further lower the lattice potential energy in order to preserve 
the mechanical equilibrium. 

to thermal phonons is given by: Ft = 
—kT\n^'^ex.p{—huJn{n-\-^)/kT), where n la- 
bels the eigenmodes. If kT ^ J^max^ one has: 
Ft = -{?>N/V)kT\n{kT/mjj), where uj is defined such 
that IncJ is equal to the average value of Incj. The 
contribution of the thermal phonons can be written as: 
Ft ^ —{3N/V)kT0j'^ where the non-dimensional factor 
= —{d'^/dj'^)^^o\nJiIU/kT has been demonstrated [9^ 
to be of order unity when the harmonic potential 
dominates the pair interaction potential. This gives a 
good estimate: Ft ~ -{3N/V)kT-f^ . 

The number of mechanical bonds per atom z requires a 
careful definition in amorphous thermal systems. If there 
were only covalent bonds, then z is obviously just equal to 
the number of covalent bonds per atom. However, in ad- 
dition to covalent bonds, weaker interactions are present, 
such as those between two monomers of different chains in 
polymer systems. Such interactions are of van der Waals 
nature and they can be modeled by the Lennard- Jones 
(LJ) potential. In that case, it is important to distin- 
guish quantitatively between these contributions in the 
total z. Here we propose the following criterion that al- 
lows us to calculate these terms unambiguously. With a 
glassy polymer in mind (although without any loss of gen- 
erality), we shall assume that a contribution to z arises 
from inter-chain interactions whenever the two monomers 
are at a mutual separation r ^ Tmin, where rmin is the 
minimum of the LJ potential well, see Fig. [2] Although 
this procedure might slightly underestimate the actual 
absolute value of the underestimation has no effect on 
the change Sz upon varying T which is the quantity that 
matters in our calculation. Later in the text, when we 
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FIG. 2: Schematic of the criterion used to define the contri- 
bution of the weaker Lennard- Jones interactions to the total 
number of mechanical bonds z. Only pairs of particles that 
lie within the soft repulsive part of the potential (to the left 
of the minimum) contribute to the z counting. 



compare the model predictions with a particular experi- 
ment on polymer glasses [20] , we shall derive the explicit 
values for the two components contributing to the total 
average z\ the one due to intra-chain covalent bonds z^^ 
and the one due to inter-chain LJ bonds, which we denote 

Having defined the total z^ we now need to relate it 
to T. Upon introducing the thermal expansion coeffi- 
cient, ttT = y {dV/dT)^ and replacing the volume V 
via (j) = vN/V^ after integration we obtain ln(l/^) = 
aTT + const (later we shall need to estimate this con- 
stant, obtaining C ~ 0.48). Now z can be estimated 
as a function of (j) by introducing the radial distribution 
function (rdf) g{r). Since the average connectivity due 
to covalent bonds remains fixed, only the weaker contact 
bonds contributing to zlj are changing upon increasing 
the packing fraction (j) by 5(j). The increment Sz can be 
calculated in full analogy with soft-sphere systems where 
only the repulsive part of the potential is active. This 
increment is given by |22J: 6z ^ il^^^ ^'^g{C)d£,^ where 
^ is the normalized center- to- center distance. The rdf 
in the repulsive part is dominated by ^(^) ^ — 1, as 
shown in theory [22 and simulations [23] . To keep things 
analytical, here we neglect the thermal broadening of the 
rdf (which could be calculated only using involved replica 
techniques [24 ). This particular simplifying approxima- 
tion should still work for relatively low- Tc polymer glassy 
systems, as we are going to verify below, but certainly has 
to be adjusted when dealing with inorganic glasses which 
have an order of magnitude higher Tc. The increment 8z 
has to be measured from the point where the system is 
marginally stable, i.e. z = Zc cii (j) — (j)ci and from the 
integral we obtain: z — Zc ^ — ^c- In the affine ap- 
proximation, the solid becomes marginally stable only in 
the limit ^ and — > 0, and hence we have z ^ 0^/^. 
Using the earlier relation between cj) and we obtain: 
z ~ g-aTT/2 Substituting z and (j) in Fa + Ft we now 
can write the full expression for the shear modulus in the 
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affine approximation, Ga = d'^{FA + -Pr)/t^7^, yielding: 

The Born criterion of melting [2] , is given by Eq. ([T]) set to 
zero: f^R^ = A:Te"^^/^ We shah see later that arT <C 1 
and, remarkably, this relation reproduces the Lindemann 
criterion [25] , which uses equipartition to state that melt- 
ing occurs when the average vibrational energy of a bond 
equals kT. The Lindemann criterion grossly overesti- 
mates melting temperatures for amorphous solids pT] : 
also, Eq.Q cannot capture the vanishing of rigidity as 
seen in the melting of glassy polymers [20 . It turns out 
that to describe the melting of amorphous solids one has 
to account for nonaffine deformations in the lattice dy- 
namics. 

The shear modulus accounting for nonaffine deforma- 
tions is derived in [26 as a lattice sum: G = Ga — Gna = 
Ga - EilXjS'j'lf where R^^ = (d'F/dddr^^o 
is the dynamical matrix of the solid. The vector /. mea- 
sures the increment of force acting on an atom in response 
to the deformation (here 7 = r]xy) of its environment. 
It can be shown that for harmonic pair potential [26] : 
/. = —Riji^^je^jefje^j^ where e^j is the unit vector 
along the bond connecting two atoms i and j. Since 
the sum runs over bonds to the nearest-neighbors j of 
the atom i, it is evident that in a perfect crystal for each 
bond involving i there is a mirror-image bond across a 
reflection plane of the crystalline lattice. Therefore, ev- 
ery bond in the sum cancels with its mirror- image, and 
/. =0, Vz in most crystal lattices. As a result, /. is 
nonzero only with lattice disorder, see Fig. ^p. The non- 
affine correction to the elastic free energy arises to en- 
sure that mechanical equilibrium, which disorder tends 
to compromise, is preserved upon deformation. It has 
been evaluated, in mean-field approximation, for random 
assemblies of harmonically bonded particles [271 IM] • 

G = Ga-Gna = ^^H^- Zc). (2) 

OTT Ho 

It should be noted that this expression has been derived 
for harmonic pair potentials and is therefore valid for any 
potentials with an attractive minimum. 

The nonaffine contribution is encoded in Eq. ([2| in 
the term proportional to Zc which expresses the internal 
energy required to fuel the nonaffine motions necessary 
for the preservation of mechanical equilibrium against the 
effect of disorder. With purely central-force interactions 
in d dimensions, the shear modulus vanishes at Zc = 2d 
{zc = 6 in 3D) because the nonaffine term is proportional 
to the number of degrees of freedom that can be involved 
in the nonaffine energy relaxation. This is consistent with 
the classical Maxwell criterion for marginal stability with 
purely central forces: G ^ {z — 6). In general, Zc defines 
the critical coordination at which the lattice is no longer 



rigid because all the lattice potential energy is "spent" 
on sustaining the nonaffine motions and no energy is left 
to support the elastic response to deformation. Using 
again: ln(l/^) = a^T + C, we arrive at: ln(^c/0) = 
arir — Tc). The corresponding relation z — Zc^ ^/cj) — (j)c 
can be manipulated into: 

\n{<l>,/^) = -\n[l + {z-Ze?/^e]. (3) 

Combining this with the relation for (pi^ar) we obtain 
ln[l + — Zc)'^ /(j)c\ = aT{T — Tc)^ and finally arrive at 
the condition: z-Zc = \/^c[e'^^'^^""^) - 1]. Substituting 
it in Eq.([2|, we obtain: 

Ga - G^A = f ^0ce"-(^-^V'^^[e""^'^^"^^-l]- 

OTT ilo * 

(4) 

According to this equation for the shear modulus G{T), 
nonaffinity alone (induced by disorder) causes the melt- 
ing at a critical point Tc with the scaling ~ V^T^^^T, even 
without the effects of thermal vibrations on the rigidity. 
Including the effect of thermal phonons in the same way 
as was done in Eq.Q, the full expression for G(T) be- 
comes: 

(5) 

The square-root cusp singularity G ~ V^^^^T predicted 
by our theory for the shear modulus has been reported 
in the most recent numerical simulations for the melting 
of colloidal glasses [29 which thus confirm the validity 
of our approach. To assess the interplay and relative 
magnitude of nonaffinity and thermal phonons, it is im- 
portant to have an expression for the critical point Tc in 
terms of parameters used in this analysis. Since Tc is, in 
effect, the glass transition temperature, the task of find- 
ing it explicitly for a general system remains challenging. 
However, since our comparison will be with a particular 
experimental system of polymeric glass, we can in fact 
offer an expression in such a case. For polymer chains 
of n units, the average connectivity due to intra-chain 
covalent bonds is close to 2: = 2(1 — 1/n) (the to- 
tal coordination number is z = z^o + ^lj)- The LJ are 
central forces, but the covalent bonds also put a con- 
straint on the bond angle. The classical Phillips-Thorpe 
analysis of marginal stability [SO] [31] gives the fraction 
of floppy modes / = A^floppy/3A/' in a purely covalent 
network: / = 1 — ^([^Zco + [2^co ~ 3]), where every z^o- 
coordinated monomer contributes 2zco — 3 bending con- 
straints, in addition to \zco stretching constraints. In 
our case, upon adding the LJ inter-chain bonds into the 
counting, the fraction of floppy modes becomes: 

/ = 1 - I {\z,o + [2^co - 3] + ^ZLj) . (6) 

Keeping z^o fixed, since it is T-independent and fixed 
by the polymer chemistry, we set / = in the above 
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FIG. 3: Comparison of theoretical predictions from Eq.([5]) 
(nonaffine + phonons), Eq. Q (nonaffine athermal), and Eq. 
|l} (affine + phonons). We take Tc = 383K, Rq = 0.3 nm, 
ar — 2-10~^K~^, as for polystyrene, (j)c — (jfc — A-Zco = 0.44, 
and the spring constant k = 50 N/m. 



equation and solve for the critical value of zlj at the 
rigidity transition, obtaining: z^^ = 12 — bz^o- Upon 
applying Zc = z^o + ^lJ' obtain the critical value of 
the total connectivity z at which the rigidity is lost: Zc = 
12 — Azco = 12 — 8(1 — 1/n) . In other words, in order for 
the amorphous polymer assembly to become solid (glass) 
there need to be at least z£j = — ^^co = 12 — 5zco LJ 
inter-chain bonds per monomer, in addition to the chain 
connectivity. Now we convert Zc into the critical volume 
fraction 0c , via (j)c = 0J — A • 2:co, where 0J is the packing 
fraction in the limit = 0. If the attraction is weak or 
absent, as in a system of hard spheres, then 0* :^ 0.64 as 
for random packings. 

Finally, using the expression for (j)c{oLTT)^ and for 
^00(^)7 we obtain the glass transition temperature for 
chains with degree of polymerization n: 



Tc = —(1 



C 



-2A) 



2A 



(7) 



The first term is what remains for very long chains 
(n ^ 2K/aT)' The expression ^ provides a theoretical 
foundation for the empirical dependence of the glass tran- 
sition on n, first discussed by Flory [32^,^33]. For common 
polymers, the experimental values of the factor 2K/aT 
[32l|33] are of order lO^i^, and hence A 0.1. If we take 
Tc — 383K, as for polystyrene glass used in the experi- 
mental comparison below, this gives a reasonable value 
of C :^ 0.48. This value implies that (/) 0.61 • e""^^, in 
dense amorphous polystyrene. 

In Fig. [3] we have plotted predictions of different the- 
ories for the shear modulus G{T) using the parameters 
of amorphous polystyrene, taking n = 200 [20j. It is 
evident that nonaffinity is the main effect driving the 
melting transition, whereas thermal phonons have prac- 
tically no effect on the qualitative behavior of the melting 
curve. On the other hand, if nonaffinity is neglected, the 
transition is shifted to unrealistic, enormously high tem- 
peratures. 
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FIG. 4: Theoretical fit of the experimental data for amor- 
phous polystyrene using Eq.([5| and the same values for pa- 
rameters as in Fig [3] including Tc = 383K, and fitting for 
the single parameter: the spring constant k — 52 N/m. The 
experimental data were taken at a constant frequency of 0.9 
Hz. 



Finally, Fig. |4] shows the comparison of our theory pre- 
diction for the melting of amorphous polystyrene with the 
classical experimental data of Schmieder and Wolf [20] . 
reported also in Ferry's monograph as representative of 
the quasistatic, low- frequency G{T) in the glassy solid 
state [34]. As discussed by Ferry, the catastrophic drop 
of G at the critical point is the hallmark of the low- 
frequency (static) response, whereas at higher frequency 
the drop becomes more gradual (due to the reduced non- 
affinity, from our point of view). Experimental data 
taken at even lower frequency would not differ much as 
the drop can hardly be sharper than in Fig. 4. 

The value of the fitting parameter, the spring constant 
hz ^ 52 N/m is very sensible: it corresponds to the C- 
C covalent bond enthalpy 350kJ/mol over a distance of 
0.15 nm, which is almost exactly the C-C bond length 
(0.146 nm). More importantly, the theory can reproduce 
the qualitative behavior of the experimental curve and 
the criticality (that are found also with many other poly- 
mers pQl[34]) very well indeed. It is important to notice 
that no alternative theory is available in the literature 
for the mechanical response of amorphous polymers at 
T < Tg (in contrast, at T > T^, the reptation theory of 
de Gennes, Doi and Edwards [35] provides a good under- 
standing of viscoelasticity of polymer melts). 

From the point of view of applications, our theory 
can be used to reconstruct the thermoelastic behavior of 
glassy polymers at T < T^, something which has not been 
possible so far, in the same way as the Williams-Landel- 
Ferry theory [34] is used for the viscosity at T > . The 
theory can be used to model also the melting curves of 
inorganic materials, such as amorphous oxides, semicon- 
ductors and metallic glasses where the thermal phonons 
play a more important role in view of the much higher 
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